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Abstract 

The production of hard photons in hadronic collisions is studied using Soft-Collinear 
Effective Theory (SCET). This is the first application of SCET to a physical, observable 
cross section involving energetic partons in more than two directions. A factorization 
formula is derived which involves a non-trivial interplay of the angular dependence in 
the hard and soft functions, both quark and gluon jet functions, and multiple partonic 
channels. The relevant hard, jet and soft functions are computed to one loop and their 
anomalous dimensions are determined to three loops. The final resummed inclusive 
direct photon distribution is valid to next-to-next-to-leading logarithmic order (NNLL), 
one order beyond previous work. The result is improved by including non-logarithmic 
terms and photon isolation cuts through matching, and compared to Tevatron data and 
to fixed order results at the Tevatron and the LHC. The resummed cross section has a 
significantly smaller theoretical uncertainty than the next-to- leading fixed-order result, 
particularly at high transverse momentum. 



1 Introduction 



The production of hard photons in high energy collisions is one of the most fundamental 
processes to be observed at any hadronic collider. Photons which are produced from the 
underlying partonic interaction are called direct, or prompt photons. They can probe the 
structure of the proton at very small distance scales, and are therefore sensitive both to 
details of the standard model and to possible new physics scenarios. At lowest order, there are 
two partonic processes which can produce direct photons: the annihilation channel qq —> 7*7 
and the Compton channel qg —> q^y. The Compton channel is particularly important as it 
gives direct access to the gluon parton-distribution function (PDF). 

On the experimental side, the photon spectrum can be measured with great precision. 
However, it is in general not possible to distinguish whether the photons are direct, that is 
they have come from the underlying hard interaction, or if have been produced from secondary 
fragmentation, such as tt° decay. This ambiguity is lessened somewhat for very high energy 
photons, which are relatively unlikely to have been produced from fragmentation. Moreover, 
demanding a mild isolation criterion on the photon, for example, that there be less than 2 GeV 
of hadronic energy in some reasonable surrounding area, makes the high px photon spectrum 
a fairly clean probe of the underlying interaction. 

On the theoretical side, the direct photon spectrum has been approached both in fixed 
order perturbation theory and with soft gluon resummation. The cross section is known 
for both the polarized and unpolarized case at next-to-leading order (NLO) [H [21 [3], in the 
fully inclusive case including analytic integration over the real emission contribution. The 
Monte Carlo program jetphox [I] implements the NLO result numerically, as well as the 
contamination from fragmentation, and allows the user to specify an isolation criterion. 

Near the partonic threshold, where the transverse momentum px of the photon is close 
to half of the partonic center of mass energy px < vI/2, the invariant mass of the recoiling 
hadronic system becomes small and the partonic cross section involves large logarithms. These 
logarithmic terms, which arise from soft and collinear radiation, often amount to the bulk of 
hadronic cross sections. To improve predictions, threshold contributions can be resummed to 
all orders in perturbation theory. For direct photon production, this has been done to next- 
to-leading logarithmic order (NLL) [51 [7J El E] and a phenomenological comparison to data 
from E-706 and UA-6 has been made [HI [ID] • The resummation effects are important at large 
Pti and therefore must be understood to improve the precision of theoretical predictions for 
the direct photon px spectrum and related observables. 

The approach we take to resummation in direct photon production is based on the use 
of effective field theory techniques. Effective field theories are powerful tools for separating 
physics associated with different scales and resumming large logarithms of ratios of those scales 
through the renormalization group. In this paper, we apply Soft-Collinear Effective Theory 
(SCET) [TTJ [121E3] to direct photon production. The effective theory was originally developed 
to analyze -B-decays, but its promise for collider physics was envisioned early on [13]. The 
collider applications of SCET have included deep- inelastic scattering (DIS) [151 EEl CEO HE] , 
Drell-Yan [T9] [20] . Higgs production [211 [221 [23] , tt production and event shapes in e + e~ 
collisions [22 [23 [2S1 [23 E2E1 [21], and electroweak Sudakov resummation [HI El E2]. The 
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main result of much of this work has been to improve our understanding of the effective 
theory description of QCD, but in some cases it has led to qualitatively new phenomenology 
For example, in [28] the N 3 LL resummation of the thrust distribution was performed, leading 
to one of the best measurements of a s and a strong model-independent bound on the gluino 
mass [33]. However, there is still much work to be done in demonstrating the power of this 
effective field theory at hadron colliders, and the present paper is a step in that direction. 

All of the previous collider applications of SCET have involved processes with only two 
directions of large energy flow: in DIS, these are the proton and outgoing jet, in Drell-Yan, 
they are the incoming hadrons, and in e + e~ they are the outgoing jets. Formal expressions 
involving several directions of large energy flow are straightforward to write down [3H [35] , 
and the theory with multiple collinear fields has been used to perform electroweak Sudakov 
resummation of partonic amplitudes [HU E2], and to derive constraints on the structure of 
infrared singularities of gauge theory amplitudes [311 E3 EB]- However, so far it has not 
been applied directly to a physical process, and it is an important step to derive and check a 
factorization theorem for a cross section involving three directions. The simplest such process 
at a hadron collider is direct photon production, which explains much of the motivation for 
the current work. 

We begin in Section [2] with an overview of direct photon production, including a review of 
the relevant kinematics and a physical discussion of the factorization theorem. In Section [31 
the factorization theorem is derived with SCET. In order to achieve resummation at the next- 
to-next-to-logarithmic order, which is one order beyond previous results, we need the one-loop 
expressions for the hard, jet, and soft functions appearing in the factorization theorem. These 
are calculated in Section HI The relevant soft function, as defined through the factorization 
theorem, depends on radiation away from the direction of the outgoing jet, and we calculate it 
to one loop. The gluon jet function, necessary for the annihilation channel, is also calculated 
at one loop. Using renormalization group (RG) invariance of the cross section, known results 
for some of the anomalous dimensions, and Casimir scaling of the soft function, we manage to 
extract the anomalous dimensions of the relevant hard, jet and soft functions to three loops. 
After solving the relevant RG equations, we combine these ingredients together into a closed 
analytical formula for the resummed direct photon distribution. In Section [51 we show that 
the renormalization scale independence of the cross section implies a non-trivial cancellation 
among angular dependent parts associated with different scales. Section [6] discusses the scale 
choices and matching procedure. Finally, in Section [7J we evaluate our formula numerically, 
comparing to Tevatron data and making predictions for the LHC. 

2 Direct photon cross section 

In this section, we establish some notation for kinematics of direct photon production. Then 
we review the differential cross section in fixed-order QCD and discuss heuristically the fac- 
torization formula which we derive with SCET in Section [31 
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2.1 Kinematics 



Let the incoming hadron momenta be and P% and the photon momentum be p*t. We are 
interested in photon production at high px = Pt- Our results will be most accurate when p? 
is near the machine threshold limit, 



Ecm 
2 cosh y 



where -Ecm = a/ (-Pi + P 2 ) 2 is the center of mass energy of the collision and y is the photon's 
rapidity. p™ ax is the maximum px the photon can possibly have for a given y. Of course, 
the phenomenology of direct photon production is dominated by much smaller transverse 
momenta, but the factorization theorem will only have exact perturbative scale independence 
for pt ~ p™ ax , and for its derivation we will expand around this threshold. 

Near threshold, the recoiling radiation X must have P^ ~ p T ~ p™ ax , which is only 
possible if the mass of the recoiling radiation is close to zero. By momentum conservation, 

P? + P?=$ + P$, (2) 

and the threshold implies Ex ~ p™ ax S> \J P x . Near this limit, the recoiling radiation can be 
characterized as a jet of collinear particles with momentum pj accompanied by soft radiation 
with momentum P x = Pj + k^. As we will discuss in the next section, SCET provides a 
field-theoretic description of the associated collinear and soft partons and their interactions. 

At leading order, there are two channels for direct photon production: the Compton process 
qg — > qj and the annihilation process qq —> g^. In either case, let the incoming partons have 
momenta p^ = X\P^ and p^ = a^^- The hadronic and partonic Mandelstam variables are 

s = (P 1 + P 2 )\ t = (P 1 -p 1 ) 2 1 u = (P 2 - Pl )\ (3) 

and 

s = (Pi + P2) 2 = xix 2 s, t = (pi - p 7 ) 2 = xit, ii = (p 2 - Pj) 2 = x 2 u . (4) 

For direct photon production, it is conventional to work not in terms of the Mandelstam 
variables, but in terms of dimensionless ratios of them. 

v — 1 + - , w = . (5) 

s s + t 



We will also use the shorthand 



for compact notation. It follows that 



v = 1 — v (6) 



s = ± P \ i=-± P A u = - P A (7) 

wvv w v v 

1 Pt y Pt - y , Q , 

x x = e y , x 2 = - — -e y . 8 

w E CM v E cm v 
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At the hadron level, the event is characterized by two quantities, pt and y. At the parton 
level, it takes four, for example, {pt, y, X\, x 2 }, or {p T ,y,v,w}. 

To understand the thresholds, it is helpful to define the hadronic invariant mass 

M 2 X = P\ = (Pi + Pi - p 2 ) = s + 1 + u (9) 

and the partonic invariant mass 

m x — (Pi +P2- P~f) 2 = s + i + u. (10) 

The partonic invariant mass, mx, includes only the partons involved in the hard scattering 
process, while the hadronic mass, Mx, includes also the proton remnants. Note that while Mx 
is observable, mx must be integrated over in any measurable quantity. In the literature (e.g. 
in [S]), the above two quantities are sometimes denoted by S4 = M\ and S4 = m\. These 
quantities represent the mass of everything in the final state except the photon, at the hadron 
and parton levels respectively. At leading order in perturbation theory, where the partonic 
final state consists of a single parton, w = 1 and mx = exactly. 
In terms of pr, y, v and w the threshold variables read 

M 2 X = E 2 CM - 2p T E CM cosh y = E 2 M (1 - pt/Pt") (1 1) 



and 



In terms of mx,%i, x 2 and v , 



m x = — . (12) 

V w 



Ml = ™* + E 2 CM [(1 - x x )v + (1 - x 2 )v] . (13) 

Since partonic configurations are specified by four variables, surfaces of constant Mx(pr,y) 
are three dimensional. That is, there are three independent ways we can have the a small 
deviation from Mx = 0. It is natural to take the independent variations to be x±, x 2 and m\ 
because setting x± — x 2 — 1 and mx = forces Mx = exactly. The fourth variable, v, can 
be thought of as moving us along the surface of constant Mx- As Mx — > 0, vv — > Pt/Eq M . 
Then, to first order in 1 — x±, 1 — x 2 and m\, 

2 

M 2 x = m 2 x + ^ [(1 - Xl )v + (1 - x 2 )v\ + . . . . (14) 
vv 

This is equivalent to the threshold expansion in [5]. This form for Mx will be convenient for 
checking the SCET factorization theorem in Section [5j 

Note that the limit Mx — > automatically enforces that the reaction takes place at the 
threshold x\ — > 1, x 2 — > 1, where the leading partons carry almost all of the proton momentum. 
In contrast, taking mx — > does not force X\ — > 1 or x 2 — > 1. At the partonic level, the 
factorization theorem will resum logs of mx, which appear as a"ln m (l — w). It will also 
resum logs from the evolution of the parton distribution functions, of the form a™ ln m (l — x^), 
which are only relevant near the machine threshold. These two types of partonic logs must 
then be combined to resum logs of a"ln m Mx/Eqwi in the full observable differential cross 
section. 
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2.2 Differential cross sections 



Using v and w, the cross section can be written in the form [3] 



dV 

dydpq 



Pt 



ab 



pt 



PT 



<:<J 



dv 



PT 1 
3 CM 



dw [xifa/N^X!, fl)] [x 2 fb/N 2 (x2, (J>)] 



(;!' 



where the sum is over the different partonic channels. 

At leading order the mass of the final state is zero, w = 1, and 



d 2 a ab 
dwdv 



a ab (v)5(l - w) 



where 



a qg {v) 



Pt 



q -a s (n)— (I + v 2 ) v 



P 2 T 



d 2 d ab 

dwdv 



(15) 



(16) 



(17) 



Here, e q are the charges of the quarks and N c is the number of colors. 

At next-to- leading order (NLO), the partonic cross section acquires w dependence. It has 
the form (leaving the partonic indices ab implicit) 



d 2 



a 



dwdv 



a(v) < 5(1 — w) + a s (/i) 



S(l-w)h 1 {v) + 



1 



+ 



1 — w 
ln(l — w ) 



1 — w 



h 2 (v) 
h 3 (v) 



h±{v, w) 



The plus distributions indicate the singular behavior at NLO as w — > 1, that is, as the 
kinematic threshold is approached. An N n LO computation would lead to higher-order plus 

distributions, up to — — 1 _^~ w - > . Keep in mind that there is implicit, non-singular, w de- 
pendence in the PDF f a (xi,fi) as well. These functions hi can be found in [3]. The singular 
ones, hi, hi and h 3 , as well as the singular coefficients at NNLO are listed in Appendix IB1 
The result from effective theory, which we derive in Section [3j has the forno 



d 2 a 
dwdv 



w a(v) H(pt, v, /i) / dk J{m 2 x — (2E j)k, /i) S(k, /i) . 



(19) 



Here, H, J, and S are the hard, jet and soft functions, respectively and Ej is the energy 
of the jet. The functions H, J, and S are different in the two partonic channels. The hard 



1 The w prcfactor in this equation is a convention, but it follows from the s 1 dependence of the partonic 
cross section and conveniently cancels the w-dependence of x\ in Eq. (|15[) . 
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function comes from matching SCET to QCD. It depends only on v, since w = 1 at the hard 
scale fih ~ Pt- The jet function comes from integrating out collinear modes. It depends on 

2 

w through m 2 x = -Ezii=w_ i n particular, at leading order J{p 2 ) = S(p 2 ), which reproduces 
the 5(1 — w) dependence of the LO cross section. Expanding the hard, jet and soft functions 
to order a s will reproduce the other terms in the NLO cross section, up to terms which are 
regular as w — > 1 and correspond to power corrections in the effective theory. Expanding to 
order a 2 produces all the singular terms at NNLO, which is a new result. 

The scale \x in Eq. (fl9|) should be identified with the factorization scale since it determines 
where the PDFs are evaluated. Since the physical scales entering H, J, and S are quite 
different, any choice of fi will lead to large perturbative logarithms. To resum these, we will 
solve the RG equations for the three functions in Section H] and evolve each of them from a 
matching scale to the reference scale \x at which the different contributions are combined. For 
the matching scale for the the hard function, we choose \ih = pr- The choice of the matching 
scales for the jet function and soft functions is less obvious, and will be discussed in Section [6j 

The form of the SCET factorization theorem, Eq. (]19J) . can be understood from simple 
physical arguments. The recoiling radiation X in a high-p^ direct photon event is almost 
massless. That is, P x <C E x . Thus, this radiation consists of particles forming a jet, with 
momentum pj and additional soft radiation k^. Then, 

m x =Px = (Pj + kf = m 2 + 2k-pj (20) 

up to terms of order k 2 <C m 2 , where m 2 = p 2 is the mass of the jet. The precise allocation of 
the final state particles into the jet or the soft sector is not well defined, but the ambiguities 
give only power suppressed corrections to this relation. 

Since the jet is lightlike at leading order, its momentum can be written as pj ~ Ejrij, 
where rij = (1, n,j) is a lightlike vector. This is true up to power corrections, because n,j-p J ^> 
pj_ ^> nj-p 3 . Thus, the amplitude for producing a configuration with a particular value of 
m 2 x will be proportional to 

J(m 2 ) = J (m 2 x - (2Ej)(n,j-k)) . (21) 

This explains the Ej dependence in Eq. ffl9|) . Also, we see that the only component of the 
soft radiation which is relevant to threshold resummation is the rij ■ k component, that is, 
the piece backwards to the direction of the jet. This projection k = (n,j ■ k) also appears 
in the soft function of Eq. ( fl9l) . Our derivation of the factorization theorem in the next 
section will provide us with operator expressions for the jet and soft functions appearing in 
the factorization theorem. We will compute these functions to one loop in Section HI 

3 Derivation of the factorization theorem 

We will split up the derivation into two parts. First, we will summarize some results from 
SCET about operators and scaling relations. Then we will apply the effective theory to the 
process of direct photon production and derive the factorization theorem for the differential 
cross section. 
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3.1 Soft-Collinear Effective Theory 

SCET provides an expansion in the limit of large energies and small invariant masses. For 
a jet with momentum pj, we expand in e = mj/Ej. For direct photon production, there 
are three high-energy scales, the energy of the two incoming partons and the energy of the 
hadronic final state X. Correspondingly, we introduce three different sets of collinear fields 
associated with the directions of large energy flow and one set of soft fields, which mediate 
interactions among the different collinear directions. 

The components of momenta p^ of collinear fields associated with the i-th direction, 
whether quark (ipi), antiquark (ipi) or gluon (Af), scale as 

i-collinear: (n* -p ci , ni -p ci , p^ x ) ~ (1, e 2 , e) , (22) 

where the vectors are light-like reference vectors along the i-th direction. For each light- 
cone vector tt.^, we also introduce a conjugate light-cone vector, n^, such that n^-ni = 2. The 
field components of a collinear gluon scale in exactly the same way as its momentum, due 
to gauge invariance. Two components of collinear Dirac spinors are suppressed and can be 
integrated out, after which the spinors fulfill the constraint fiiXi = 0. Then, collinear fermion 
fields scale as e. All components of soft momenta, with respect to any of the jet directions, 
are small 

soft : p» ~ e 2 . (23) 

Thus, soft fields can interact with any of the collinear fields without changing their scaling. 
Soft gluon and soft quark fields scale as ~ e 2 and ip s ~ e 3 , respectively. 

To construct operators in the effective theory, it is convenient to work with the jet fields 
Xi and A^ ± [121 [39] . They describe the propagation of energetic partons in the i-th direction, 
but do not correspond to any experimental definition of jet, such as cone or /cr-jet. Explicitly, 
the jet fields are related to free collinear quark and gluon fields by the addition of Wilson lines 

Xi (x) = W}{x)^^{x) , A1(x) = W}{x) [iD$Wi(x)] . (24) 
These z-collinear Wilson lines 

Wi(x) = P exp ^ig J dsni-Ai{x + snj)J (25) 

ensure that fields are invariant under collinear gauge transformations in each sector [TT| fl~2]. 
The symbol P indicates path ordering, and the conjugate Wilson line W} is defined with the 
opposite ordering prescription. 

At leading power only the rii-A s component of the soft field can interact with the collinear 
fields in the i-th direction, since all other components are power suppressed compared to 
components of the collinear gluon field. As a consequence, the leading-power soft- collinear 
interactions are Eikonal and the soft dynamics can be removed from the collinear Lagrangian 
to all orders in perturbation theory through field redefinitions. For example, the interaction 
of soft gluons with collinear fermions in the SCET Lagrangian has the form 

C Ci+s = Xi(x)^n i -A s (xJ)xi(x). (26) 
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where xt = (fii-x)^-. The peculiar x-dependence of the soft field will be explained below. 
This interaction can be represented in terms of soft Wilson lines. Redefining the quark and 
gluon fields as 



where 



Xi{x) -> Yi(xJ) Xi(x) , 
Xiix) -»■ Y?(x-)xi(x) , 
^ ± (*)^(*_Mf ± (*)r/(*_), 

Yi(x) = P exp (ig [ dtn, r A a s (x + tni)t a 



(27) 
(28) 
(29) 

(30) 



eliminates the interaction C Ci+s and other pure-gluon terms. After this decoupling transfor- 
mation, soft interactions manifest themselves only through Wilson lines in the operators [12] . 

Let us now explain why the collinear fields in soft-collinear interactions, such as in Eq. f )26|) . 
are evaluated at x and the soft fields at x_. First, recall that a collinear sector alone should 
be equivalent to full QCD (it can be derived as QCD in a boosted frame). Therefore, no 
information must be lost in the derivative (or "multipole") expansion of a collinear field 



ipi(x) 



1 + ~{rii-x)d ni + -(rii-x)d ni + x±-d A 



^i(O) = i>i{x) + 0{e) 



(31) 



(e 2 , l,s), for each of these terms to to scale like e° the scaling of x is 



Since (d nv d nv d A 
fixed: 

(ni-x,ni-x,x±) ~ (e~ 2 , l,^" 1 ) . (32) 

This is collinear scaling in position space. Now, for a soft field interacting with a collinear 
field at the same point x, we can multipole expand as well 



1 l 

1 + -{n i -x)d ni + -(ni-x)dn t +x±-d± + 



A!(p)=AZ[(n i -z) 



n- 



0{e) 



(33) 

For these fields, since the soft momenta scale like (d ni , d Hi , d±) ~ (e 2 ,e 2 ,e 2 ), only the terms 
like \{fbi ■x)d nj ] k are unsuppressed, which is why A%(x) = A^(xJ) at leading power. This is 



simply the position space version of Eq. (|20|) . (pj + k)' 



m 



2Ej(nj-k) + 0(6), and will 



play a crucial role in the derivation of the factorization theorem below. 



3.2 SCET for direct photon production 

To study direct photon production in the effective theory we first introduce three light-like 
reference vectors. Two vectors and are aligned with the beam and point in the direction 
of the incoming hadrons with momenta and P%. The third reference vector rij is along the 
direction of the hadronic jet, which recoils against the hard photon. 

With the field content and scaling dimensions of SCET established, the first step is to match 
to the full standard model. For direct photon production, we need to introduce operators which 



S 



can reproduce the matrix element of the the vector current J u (x) = ip(x) ^ v ip(x) in the full 
theory. The leading operators relevant for the partonic process qq — > 7$ are 

O s q y(x»; h, t 2 , tj) = x 2 (x" + t 2 n%) Aj v ± (x» + tjn*) X i(x" + tin?) , (34) 

T q y{x»- h,t 2 , tj) = x 2 (^ + t 2 n^) %a vp A/ ± (x» + f jf#) X x{^ + ti<) . 

These two operators correspond to the two cases where the quarks have equal or opposite 
spin. We are interested in the unpolarized cross- sect ion, and so will need the sum of both 
contributions. In addition to O qq ? and there are 5x2 more operators, which are obtained 
by permuting the indices on the fields in O qq u and O^, and which describe the other partonic 
channels with initial states qq, qg, gq, qg, and gq. Their Wilson coefficients can all be derived 
from the Wilson coefficients of the O qq operators by exchanging the momenta. For the case of 
the O qg operators, associated with the qg — > jq channel, the corresponding crossing relations 
are nontrivial, and care has to be taken to get the proper imaginary parts. There are also 
operators with three collinear gluon fields. It is straightforward to include them, however, the 
99 ~^ 91 channel starts contributing only at NNLO and these operators are only relevant for 
N 3 LL resummation. In the following, we will generically refer to all the operators relevant for 
direct photon production as CJ. 

Note that the operators Oj are not local. The non-locality arises because derivatives 
along the directions associated with large momentum flow are not suppressed. The variables 
tx,t 2 and tj on which the operators depend are the position space equivalent of the label 
momenta introduced in [TT]. The smearing in the direction which they induce allows for 
different amounts of energy in the corresponding collinear fields. The Wilson coefficients for 
the operators must also also depend on these ti's, and these variables must be integrated over 
in matching to the full theory 

r{x*) =Y,f tiito 2 dtjC j (t 1 ,t 2 ,tj)OZ(x' t ,t 1 ,t 2 ,tj). (35) 
3 

For the calculation of the cross section, we will need the Fourier transformed Wilson coefficients 

C(n 1 -P u n 2 -P 2 ,nj-P x ) = J dtidtadtje-^^^^^^'M^i,^,^) (36) 

which depend on the large component of the momenta in each of the three directions. The 
fact that the Wilson coefficients depend on the large light-cone components of the collinear 
particles is characteristic for SCET. An alternative to the position space formalism [13] we are 
using is the label formalism [TT], where the large momentum component is treated as a label 
on the collinear fields, similar to the heavy quark velocity in HQET. 

The starting point for the factorization theorem is a generic expression for the cross section 
in terms of matrix elements of the production current, summed over final hadronic states, and 
differential in the photon momentum 

^ = p2 e q (9 J; F ^vr) 4 ^ + P * ~ Px-Py)\(* I e»J"(0) \Ni(Pi) N 2 (P 2 ))\ 2 . (37) 

^CM \ Z7V ) Z£j l x 
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The states \X) are the hadronic part of the final states allowed in the process and e v is the 
photon polarization vector. For high-p T direct photon production, these states must include 
a hard jet and so the scaling of P^ is like that of the jet momentum pj. 

For the matching step, we integrate out the hard modes of the theory. This amounts to 
plugging in the representation of the electromagnetic current operator in the effective theory 
and restricting the final states to soft and collinear modes. After matching the current J u (0) 
using Eq. (13"5]h the collinear fields in Eq. (137]) are evaluated at positions XiiU^D- These points 
can be translated to ti — 0, using 

X i{m$) = e +i ^'^ Xi{V)e~ iti{flvPi) • (38) 

Then we can have the momentum operators act on the states, where they evaluate to the large 
momentum associated with the given direction. Performing the integral over the convolution 
variables ti, as in Eq. (|36|) . yields the Fourier transforms of the hard matching coefficient. 
Thus, we have 



2>Ka e e 2 q d 3 p 7 \- f r]4xei{Pl+ p 2 _ Px _ Pi)x 



X 



J2C 1 (n 1 -P 1 ,n 2 -P 2 ,n J -P x )(X\e u O^0)\N 1 (P 1 )N 2 (P 2 )) , (39) 



j 



where Oj(x) = 0, 0, 0). Note that this has been separated into a sum over states \X) 

with only soft and collinear fields, and a sum over the operators. 

We would like to get rid of the sum over states and write the above expression ( 139]) as a 
forward matrix element. To this end, we turn the exp(iPjx) factors to operators exp(iPjx), by 
using the states |A^(Pj)) and \X). Then these operators act on the the three collinear fields in 
Oj(0), moving the entire operator to Oj(x). After summing over photon polarizations, this 
gives 



= F / q (0 * J2 E ^ifa ■ P ^ ■ P2 ' • ^CKn, ■ P u n 2 ■ P 2 , nj ■ P x ) 

x Jd 4 xe- l ^(N 1 (P 1 )N 2 (P 2 )\O;\x)\X)(X\O^0)\N 1 (P 1 )N 2 (P 2 )). (40) 

This has the form of a two point function, where the interaction between 0^\x) and O%(0) is 
mediated by exchange of final state particles \X). Since P x scales like a collinear field in the 
rij direction, x must scale like the conjugate position space coordinate (see Eq. (1321): 

(nj-x,nj-x,x±) ~ (e~ 2 , . (41) 

This scaling will help define the soft function below. 
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Next, we perform the field redefinition to remove the soft interactions from the Lagrangian. 
The operators then factorize into a soft and a collinear part. For example, 



° S J = (ft A T±Xi) (YlYj^YjY,) = O c q »O s qq , (42) 

where we suppress the color indices of the quark fields, which are contracted with those of 
the soft Wilson lines. The collinear operators Oy have the same form as the operators in 
Eq. (|34p . but are built from fields which no longer have soft interactions. The matrix elements 
of the collinear operators give the PDFs and the jet functions. They are sensitive to the 
gluon's polarization and quark spins, but diagonal in color space. On the other hand, the soft 
interactions are independent of spin, but inherit their color from the original process in full 
QCD. Explicitly, the soft part of the operators for the qq — > jg process are given by 

[O s Jx)] a = \Y^x)Y J (x)t a Y}(x)Y 1 (x)} , (43) 
J *J L Uj 

where i and j are color indices. Each of these Y's is a matrix in color space, and the final soft 
operator depends on the color of the quarks and the gluon. The fact that the collinear matrix 
elements are color diagonal implies that in the matrix element squared the color indices of 
Oj\x) get contracted with C^(0) so that the soft function will involve a color trace. 

Because the soft and the collinear sectors no longer interact among each other, the matrix 
elements of the operators factorize into a product of matrix elements. Also, since P% scales 
like pj, the states \X) have collinear radiation in the rij direction and soft radiation, but not 
collinear radiation in the direction of the nucleons. For the matrix element of O qq v this means 



{N 1 {P 1 )N 2 {P 2 ) O s j\x)O s qq \0) N 1 (P 1 )N 2 (P 2 ) 

(N 1 (P 2 )\xi a (x)xw(0)\N 1 (P 1 )) x (iV 2 (P 2 )|^(x) X 2a(0)|iV 2 (P 2 )) 

x ^(0|^ J± (x)|X c )(X c |^ J± (0)|0) x J](0|<-(x)|X s >(X s |C»^(0)|0> . (44) 

Note that the Dirac indices a and /3 are contracted among the different collinear fermions. 
In this factorized form, it is now obvious that the collinear matrix elements are diagonal in 
color space. As stated above, this implies that the color indices of the soft operator shown in 
Eq. (14"3"]1 are contracted between O s qq {x) and O qq (0). 

The matrix elements of the collinear fields associated with the jets give rise to the quark 
and gluon jet functions 



xM*)r x S(o)|o> = ^tr 



Ar ± (^ b jl(0)\0} = 6 ab (-gT) 9i 



d 4 p 
(2^)3 

d 4 p 



(27T) 



9(p°) [n J -p)J q {p 2 )e- i ^ 
;0(p°)J g (p 2 )e- lx \ 



(45) 
(46) 
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where i,j and a, b are the color indices of the fields. The quark jet function is the same 
universal object appearing in the factorization theorems for DIS [T7] and thrust [2j)J [2JJ 125] . 
The gluon jet function has appeared in the analysis of quarkonium production (101 EIJ H2]- 

For the other matrix elements, we can simplify things by using the SCET scaling relations 
to project out the leading power. For the soft operator, using the multipole expansion and the 
scaling of x we find that it depends only on x^L = (nj-x)^j- to leading power, as in Eq. f l33j) . 
Then the soft function relevant, for example, for the partonic channel qq — > jg reads 



(YjYjt-YfrM |0> 



dk + e-'M^v-Vs S q g(k + ) . 



(47) 

Since the soft function only depends on x_, its Fourier transform only depends on k + = rij ■ k, 
where /c M is the soft radiation in \X). The soft function can also be written as 



1 



Y^0)Yj(0)t a Yj(0)Y 2 (0) 







(2n)5(nj-p Xa - fc+) 



(48) 



where the color indices of the Wilson lines need to be contracted as in Eq. ( I4"7j) above. The 
soft function is the amplitude squared for the emission of a set of soft partons from the three 
Wilson lines. The time-ordered product appears because cross sections are extracted from 
expectation values of time-ordered products of fields. A discussion of how expressions such 
as Eq. f)47p arise in the path integral formulation of SCET can be found in Appendix C of 
Ref. PDJ- 

For the matrix elements involving the incoming nucleons, the momenta and derivatives 
scale like 

(d ni , d^, d ± ) ~ (m ■ Pi, n t ■ P i7 P ± ) ~ (e'\ 1, e') , (49) 

where e' = / Ei. We assume that the nucleon masses are negligible, mj, so that e' 

e. Then, including only the leading power in the multipole expansion, Xi{ x ) = X% ^Var^M .@ 
The expanded collinear matrix elements are the usual PDFs 



(Ni(Pi)\xi I 



-u 

K 



r Xl (0) \Ni(Pj) = -r^Pitr^r] / d£f q/Nt (Z)e 



i€(ni-x)(n i -P i )/2 



rij ■ x- 



i C 



fg/NAO e 



i£(ni-x) (ni-Pi)/2 



(50) 



for quarks and gluons, respectively. The SCET matrix elements are identical to the PDFs 
defined in QCD because the collinear Lagrangian is equivalent to the original QCD Lagrangian 
after the decoupling. Negative values of £ correspond to the anti-particle PDF, fq/Ni(C) — 



2 A proper treatment of a theory with two expansion parameters £ and e' would involve messenger modes, 
i.e. soft modes involving the expansion parameter e'. In this case, the messenger modes can be absorbed into 
the parton distribution functions. A detailed analysis of an an analogous situation has been performed for 
DIS in [17], and we choose not to repeat it here. 
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fq/uM) = -fq/Ni{-£) and f g / Ni (C) = fg/uM) = -fg/Ni{-£)- Matrix elements which involve 
different fields all vanish. For example, 



(0\Xj(x) A^(0)\0)=0. (51) 



Furthermore, because of the traces in the collinear matrix elements, Eqs. (1451) and ( 15011 . the 
mixed tensor-scalar matrix elements vanish as well. Thus, only the diagonal terms j = k 
contribute in the sum, Eq. (1401) . 



Now let us combine the different ingredients. For the operator, we get 



dacxy d 4 x J deijf d6 j dfc + |cji 2 / gM (ei)/^ 2 (e2)j 9 (^)^-(A: + ) 

x e -i(P~ l x) e -ik + (nj-x)/2 e -i(p J -x) e i£ 1 (n 1 -x)(n 1 -P 1 )/2 e i&(n2-x) (n 2 -P 2 )/2 

The x integral gives {2-k) a 8^\pi — P2 ~ P^ ~ Pj ~ h+ I 2 ~)j where the parton momenta are 
Pi = £i(ni • -Pi) ^9" and p^ = £2(^2 • Pz)-?- Doing the p j integral then gives 

dace dxi / dx 2 / dk\C qq \ 2 f q / Nl (xi)f q/N2 (x 2 )J g (m 2 x - (2Ej)k)S qq (k) , (53) 
Jo Jo J 

where m x = (p 7 — p\ — P2) 2 , 2Ej = fij ■ {j>\ + p 2 — p 7 ), and we have relabeled as X{ and fc + 
as k. 

To get to the final form of the factorization theorem, we observe that at leading order 
J{p 2 ) = o{p 2 ) and S{k) = S(k). Thus, the sum over Wilson coefficients |C||, including the 
factors of 2 and such from the T traces and n-n factors, must reproduce the full leading order 
direct photon cross section. So, we define hard functions H qq and H qg for the two channels 
to be the the sum over the squares of the relevant Wilson coefficients normalized to their 
values to leading order in perturbation theory. Including the appropriate Jacobian factors, 
the contribution of the annihilation channel to the cross section reads 



d 2 a qq 2 f 1 Ec M e 



dv / dw[(wx 1 )f q/Nl (x 1 ,fi)] [x 2 f q /N 2 (x 2 , fi)] 



dydpr PT JJ^-eV J-p-^eV 

fc CM b CM v 

x a qq (v)H qq (p T ,v, ll) J dkJ g (m 2 x - (2Ej)k, Li)S qq (k, li) , (54) 



and for the Compton channel 



d 2 a 2 [ 1 ~^cM e ~ v f 1 

j j qg = — I dt> / dwUwx 1 )f q/Nl (x 1 ,Li)][x 2 f g/N2 (x 2 ,Li)] 

dydpr PT J^-eV J^^ev 

E CM E CM v 

xa qg (v)H qg (p T ,v,Li) J dkJ q (m x - (2Ej)k, /x)S qg (k, li) . (55) 
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The definition of a and a simpler, more physical, discussion of this factorization formula were 
given in Section I2T21 

To finish, let us briefly discuss the other photon-production mechanism, where the photon 
is produced by fragmentation. In this case, the relevant SCET operators involve four collinear 
fields, in the directions of the incoming hadrons and the outgoing jet as well as in the direction 
of the outgoing photon. The matrix element of the collinear fields in the photon direction 
corresponds to the fragmentation function. Since the invariant mass of the hadronic final 
state is small near threshold, it cannot contain any hard collinear partons in the photon 
direction. The outgoing collinear quark must thus fragment into the photon and a soft quark. 
Soft quark fields are power suppressed, which explains the smallness of the fragmentation 
contribution at large pr- 



4 Calculation of the cross section in SCET 

With the factorization formula in hand, we can proceed to calculate the hard, jet and soft 
functions in perturbation theory. Then we will use the RG to run between the relevant 
matching scales providing the final resummed distribution. 



4.1 Hard function 

The hard functions H qq - and H qg entering the factorization formulas, Eqs. fl54"|) and fl55|) . are 
given by the absolute value squared of the Wilson coefficients of operators, such as 0^ v in 
Eq. ( 13"4"|) . which are built from three collinear fields along the three directions defined by the 
beams and the outgoing hadronic jet. The Wilson coefficients of the operators are determined 
by calculating the qq — > 7(7 and qg — > 'yq amplitudes in SCET and in QCD. The matching 
calculation is greatly simplified by the fact that all of the on-shell SCET diagrams are scaleless 
and vanish in dimensional regularization. In the MS subtraction scheme, this allows us to 
directly read off the result for the Wilson coefficient from the fixed order calculation in the full 
theory. To this end, we use the paper [13] which gives the result for the virtual corrections to 
qq — > '-/g and qg — > 79 at one loop. In the effective theory, this result corresponds to the bare 
Wilson coefficient squared. After normalizing to the tree-level and removing divergences by 
renormalization, we then obtain the result for the one-loop hard function. For the annihilation 
channel, the result is 



Air 



2 2 
(2CV + C A ) In 2 ^ + (4CV Hvv) + 6C F ) In ^ 



+ 



-336 + 65tt 2 17 



18 



hit; lnu H — ln 2 (uu) 
6 



11 



\rt(vv) 



3 + 2w)ln 2 t;+ (48^-26) In v + (22 - A8v) \nv + (-1 -2v) InV 



6^ 2 



v 2 ) 



(56) 
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The Casimirs in the second and third lines have been set to Ca = 3 and Cp = 4/3 for 
simplicity. The result for the Compton channel is presented in Appendix [A] 

To perform resummation, we need the anomalous dimension of the corresponding SCET 
operator. This anomalous dimension is linear in In /i, which is characteristic for problems 
involving Sudakov double logarithms. For NNLL accuracy, we need the logarithmic part of 
the anomalous dimension to three loops and the remainder to two-loop order. The anomalous 
dimension of a general leading-power SCET operator for an ra-jet process involving mass- 
less partons was given in [361 EZ] and the result has been generalized to the massive case 
in [381 HH S3- The anomalous dimensions of SCET operators are related to infrared singu- 
larities of QCD amplitudes [36]. A two-loop formula for these divergences was proposed by 
Catani [46]. However, he did not have a result for the 1/e pieces at two loops. His formula was 
later derived in [47] and the missing piece was related to a soft anomalous dimension which was 
calculated to two loops in [4B1 EH] . Recently, it was realized that there are strong constraints 
on the infrared divergences, in particular from soft-collinear factorization and collinear limits 
of amplitudes [361 Effl E7J E] . These constraints explain the two- loop result for the soft anoma- 
lous dimension obtained earlier [481 EH] an d i n our case completely determine the anomalous 
dimension to three loops. 

We need the result for the three-jet operators of the form xi •^iiX'i^ as m i n Eq. fl3~4|) . In 
the color-space formalism [521 [53] used in these papers, the RG equation has the form 

\C{{Pq,Pq,Pg},fj) = T ({ P } , fJ,)\C '({Pq, Pq, Pg} , fi) (57) 



d ln/i 



E^7c Usp (a s ) In^f + £y 



a. 



\C({Pq, Pq , P g},fl) 



where Sy = la^pi ■ pj + iO, and the sign factor = +1 if the momenta pf and are both 
incoming or outgoing, and = — 1 otherwise. The Wilson coefficients only depend on the 
large components of the momentum, so p% — > |(fij -pj)nf, where = ni, ri2 or nj is the 
light-like reference vector in the direction of the appropriate parton. The color-generators are 
(Tq) a p = t^s, (Tq) a p = —t a p a , (Tg)bc = —if abc . The anomalous dimension coefficients entering 
the above equation were given to three loops in [37]. The single-parton terms involving 7* 
depend only on the representation of the zth parton and are given by two anomalous dimensions 
7 9 and 7 9 . Note that these anomalous dimensions are different from 7^ and 7^ 9 , which are 
relevant for the evolution of the PDFs near the end-point (see Eq. ( 1961) below). 

The above form Eq. ( 1571) is exact at least up to three-loop order. Terms involving the 
conformal ratios introduced in [50] can only appear for four or more partons and an additional 
constant term is ruled out by considering constraints from collinear limits [37]. Furthermore, 
for the operators we consider, there is only a single color structure: the three fields are 
contracted with where i,j and a are the colors of the anti-quark, quark, and gluon fields 
respectively in the operator. The Wilson coefficient in color space can thus be written in the 
form 

\C({Pq,Pq,Pg},fi) = t^C({p q , Pq, Pg} , fi) . (58) 
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Plugging in the explicit form of the generators, the RG-equation becomes 



dln/i 



C({Pq,Pq,Pg},V) = { 7cus P (a 



-iC P - 



Ca 
' 2 

C A 



In 



+ In 



In 



— s 



'/.</ 



+ 2f + f (59) 



The Wilson coefficients C depend on the directions as well as on the momenta pf. 
However, the dependence only arises via the large momentum components, — > ' Pi) n i ■ 
At leading power, products of these large components are equal to the usual Mandelstam 
invariants. That the hard function only depends on these invariants is also clear since it 
arises from a calculation entirely within the full theory, which has no access to the light-cone 
reference vectors, so that we know that the final answer can only depend on Lorentz-invariant 
products of the momenta. Moreover, there is only one dimensionless ratio at threshold, so we 
know that H can only depend on v = 1 + t/s. From the above result for the RG equation of 
the Wilson coefficient, we then find that the hard function for the qq satisfies 



dH q q(p T ,v,n) 
dln/i 



2 

(2C F + CU) 7cusp ln^- 2CV 7cusp ln(w) + 2 7 H 



a. 



Hqq(pT,V,n) , 



(60) 

where 7# = 2 7? + j g . The extra /3(a s ) piece comes from the our normalization of the hard 
function; it compensates for the scale dependence of the a s factor in the leading order cross- 
section (see Eq. (TIT]) ). The solution is 



H qq {p T ,v,n) 



a s (fi) 



exp [(4C F + 2C A ) S(fx h , fx) - 2A H (fx h , ft)} 



x 



-(2C F +C A )A r (^, M ) 



(vv) 2CFA ^^H g g{p T ,v,fi h ), (61) 



where H(pr, v, fit) has the perturbative expansion in a s given in Appendix El For the Comp- 
ton channel, 



Hqg(pT,V,fi) 



a s (fi) 



exp [(AC F + 2C A ) S(ft h , fi) - 2A H (fi h , fi)\ 



p% 
1% 



-(2C F +C A )A r ^ h ,fi) 



Hqg{p T ,V,fl h ) 



(62) 



The functions S(u, fi) and A(v, fi) are the same as in previous papers [T7J [28] , with a factor 
of Cp factored out of the cusp anomalous dimension in S(u,fi) and Ax{v, fi). That is 

A T (u, fi) 



S(v,fJ>) = - 



7cus P (a) 



^ , 7cusp(«) 
a B (u) P[ a ) 



(63) 



L a {u) P(a') 

Ah(v, fi) is the same as Ay but with ^ H replacing 7 CU s P - Explicit expressions for these functions 
in RG-improved perturbation theory can be found in [T7] . 
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Figure 1: Diagrams contributing to the soft function at NLO. 



4.2 Soft functions 

We consider the soft functions next. The Lagrangian of the soft sector of SCET is identical 
to the standard QCD Lagrangian, so the calculation of the soft matrix element is the same as 
in QCD. They are determined by matrix elements of time-ordered products of three Wilson 
lines. Rewriting Eq. ( 1471) . for the two channels, 



— {0\TrT[(yiYjrY}Y 2 )(x- 
±r(0\TrT[{Y}Y 2 rYlYj){x. 



(Y 2 Yjt a YjYx)(0) |0> 



(XjY^Y^O) |0> 



dk + e -^+(n,-x)/2 Sq g(k + ) 



ik + (nj-x)/2 g 



The soft functions for the qq and qg channels differ only by which representation of SU(3) 
is associated with which direction. In particular, the position xt = (fij-x)-^- at which they 
are evaluated points in the direction of the adjoint in the qq —> #7 case and a triplet (or 
anti-triplet) in the qg — > q-y case. 

In dimensional regularization the virtual graphs contributing to this soft function vanish, 
so we are left with real emission diagrams. These can be drawn as cuts through diagrams 
with a gluon being exchanged between any Wilson line at and any other Wilson line at x, 
as shown in Figure [TJ The soft (Eikonal) Feynman rules give a factor of , for the emission 
from leg i, so in particular graphs involving emission and absorption into the same leg vanish. 
As indicated by the one- dimensional Fourier transforms in Eq. (1641) . the x_ dependence means 
we only need the dependence on the component of soft radiation backward to the direction of 
the jet. 

The non- vanishing diagrams for the qq — > g~f case give 



S qq -(k) = 2 <? s y £ J j^ r - 1 6(q 2 )e(q )5(k -nj-q) 



2 J (ni-q)(n 2 -q) 



and for the qg — > q^ channel 



1 g rij- ni 

2 {n,j-q){ni-q) 



\ c rij- n 2 
2 (nj-q)(n 2 -q) 



(64) 
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S qg (k) = 2 g]^ J j0z-AQ 2 nq O )S(k -nj-q) 



1(7 ni-n 2 
2 {ni-q)(n 2 -q) 



Cp C A 

2 



nj-nx 



(n J -q)(n 1 -q) 

So the calculation boils down to the evaluation of the integral 

n a -n b 



1(7 n J' n 2 
2 {n,j-q){n 2 -q) 



I s (k)=^ J d d qS(q 2 )e(q 



[n a -q)(n b -q) 



5(k -n c -q) 



(65) 



(66) 



which we need both in the case when n£f = and in the case when is different from both 
and n%. 



To evaluate the integral, we write 



n a b n ab 



with n a ■ q± = n b ■ q±_ = and riy = n» ■ rij. In this basis, 



n 



—K H ni + rr c± 



n a b 



(67) 



(6? 



Rewriting the phase-space integration as an integral over the light-cone components and inte- 
grating over \q±\ in d = 4 — 2e dimensions, we find 



T (U\ _ 1 ,.2s ( n ac n bc 



ah 



dQ 



d-2 



dg + dg 



(l+Q- 



q~ - q+ + 2y/q + q_ cos( 



(69) 



where 9 is the angle between nj: and q±. The prefactor shows that unless the three light- 
cone vectors are distinct the integral is scaleless and vanishes. Parameterizing q + = kyx and 
g_ = kyx = ky(l — x) and integrating over y gives 



2n ab 



2s 



f,l+2s 



dxx 



1 + 2a/xxcosi 



\2e 



(70) 



We can use that the integral is symmetric under x — > x to integrate only from x = . . . \ . 
Divergences then appear only in the integration region around x = 0. After rescaling x — > x/2, 
the integral can be expanded in e, using the fact that 



x 



-l-e 



: 8(x) + 



"1" 




lnx 




— e 




X 


+ 


X 



The result is 



2tt //i\ 2e / n ac n 6c 



k \k. 

To expand this result in e, we use 



i(f)^->- 



"1" 


[fc,fi] 




"1 




[M 






-2e 








A 






A- 


//_ 


* 



+ 



(71) 



(72) 



(73) 



These star-distributions are like plus distributions for a dimensionful variable [21] • The su- 
perscript, as explained in J2H], makes explicit the singular variable and the upper limit of 
integration: these distributions vanish when k is integrated integrated from to fi. 
Putting everything together, the soft functions S qq - and S qg to order a s are 



Si(k,ii) = 6{k) + 



An 



a 



2 In 2 



2n 



12 



7T 



nijn 2J 



S{k) + 



16 In 



2n 



12 



fl y nijn 2 j 



(74) 



where the color factors for the two channels are 

1 

C q q = Cp — -CU and C qg 



:C A . 



(75) 



Note that the In 2 2 " 12 term would be absent if we rewrote the star-distribution entirely in 

nijri2j J 

terms of k\ 2 " 12 rather than k. 

ni.;n 2 j 



To get the higher order soft function we use the constraint of RG invariance to express the 
soft anomalous dimension in terms of the other anomalous dimensions (see Eq. (11041) below): 



rySqg _ 



r^ J 1 _|_ ryJg _j_ ry] 



(76) 



Since the three-loop hard, quark-jet, and PDF anomalous dimensions are known, this gives us 
7 5 ™ to three loops. Casimir scaling, which is known to hold up to at least three loops, then 
determines the other soft function anomalous dimension to three loops as well: 



2C F -C 



c. 



(77) 



The resummation of the soft functions, from fi s to \i follows from solving its RG equation 
in Laplace space, as in [TB]. The result is 



Si(k,fj.) = exp[-ACiS(fjL a ,fi) + 2A Si (ii s , ii)]si(d v! 



2n 



12 



r 



k \Hs V nun-ij) T(r], 



(7? 



with 



The Laplace transform s(L) is given in Appendix [Al 



4.3 Jet functions 



The quark jet function is known completely to two loops [55], and its anomalous dimension 
to three loops [17]. To order a s , it is 



Up\ri = S( P 2 ) + (^j {[C F (7 - tt 2 )] 5(p 2 



4C F ln^ -3C f 

p2 



(79) 
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Figure 2: Diagrams contributing to the gluon jet function at NLO. The usual gluon self-energy 
contributions are represented by the first graph. In the remaining diagrams gluons are emitted 
from one of the Wilson lines, which are denoted by crosses. 



The only place where the gluon jet function has appeared previously is in the analysis of 
quarkonium production [iQl HU H2] . In [32] , its one- loop anomalous dimension was calculated. 
Here, we will compute the full order a s gluon jet function and derive its anomalous dimension 
to order af, although for NNLL resummation we only need the a 2 s result. 
The gluon jet function is defined by 



(0) |0> 



-9?) 8*9*. 



d 4 p 

~{2TTf 



e(p°) J 9 (p 



2^ e -ipx 



(80) 



The strong coupling constant g s on the right-hand side is the bare coupling; the collinear 
gluon fields were defined in Eq. (T24|) . These collinear gluon operators only have non- vanishing 
matrix elements for intermediate collinear states. Thus, this jet function can be thought of as 
the result of integrating out the collinear modes at the scale Equivalently, we can extract 
the jet function from the imaginary part of the time-ordered product of collinear fields 



1. 



-Im 



TV 



d 4 xe ipx > 



<0|T{4£(4A£(0)} |0> 



-gT) 6 ab g 2 s Up 2 



This second definition shows that the jet function is given by the imaginary part of the 
Feynman propagator in light-cone gauge, since in this gauge the Wilson lines in Eq. (j2l|) are 
absent. 

The relevant diagrams in SCET are shown in Figure [2j In Feynman gauge all of the graphs 
in the bottom row vanish. The first graph contributes to the wavefunction renormalization. 
Since the collinear sector of SCET is equivalent to full QCD, this graph can be found in 
textbooks. In units of the tree- level result, the graph gives 



J- n 



Air 



~P 



\ Ca - \T F nA \ + ^C A - ^T F n f 



(82) 



The second and third diagrams have been computed in [56] and [57] in Feynman gauge. They 
give 



4tt 



-p z 



2 1 7T Z 

_ + - + 2- — 

e l e 6 
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Adding the contributions of both diagrams and taking the imaginary part of the propagator 
using the relation 

Im [(-p 2 - ieY l ~ £ ] = - sin(vre)(p 2 )" 1 - e , (84) 

we obtain the bare NLO gluon jet function. The jet function can be expanded in distributions 
in analogy to the soft function, see Eq. (!73|) . Renormalizing the jet function in the MS scheme, 
we then find 



C 



67 
~9~ 



7T 



20, 



-T F n f 



S(P 2 



AC A \n^- (3 
P 2 



35) 

We have also computed all graphs in gauge and in light-cone gauge and have checked 
that the result for the gluon jet function is gauge invariant. In light-cone gauge, the Wilson 
lines Wj(x) defined in Eq. ( 125]) are equal to 1 and the calculation is particularly simple, since 
only the first diagram contributes. The free light-cone propagator is 



GW(p) 
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Note that the Mandelstam-Leibrandt (ML) prescription to regulate the n-p — > singularity 
is not appropriate for SCET. The ML prescription 
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n-p 



n-p n-p + ie 
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cures the collinear singularity in the propagator, but in our case this singularity has a physical 
meaning. The Wilson line and the associated light-cone propagators arise from expanding 
QCD diagrams around the large-energy limit and the choice of the +ie prescription is dictated 
by the QCD diagrams. The loop integrals contributing to the jet function are unambiguously 
defined in dimensional regularization. They depend on a single four momentum p^. If a given 
loop integral involves m light-cone propagators it scales as (n • p)~ m so that the final result 
of the integral is independent of the sign of the ie prescription adopted for the light-cone 
propagator. 

To get the gluon anomalous dimension to higher order, we note that RG invariance for the 
qq — > jg channel implies (see Eq. f ll03f) below) 



2 7 J 



7& + 7* 



Since the PDF, hard, and soft anomalous dimensions are known to three loops, this fixes the 
gluon jet anomalous dimension to three loops as well. The result is given in Appendix lAl 

The RG equations for the jet functions from fij to [i are solved with Laplace techniques [i~6] . 
The results are 



1 fp 2 

Mp 2 ^) = exp[-4C F 5 , (/i j ,/i) + 2A Jq (ii j ,fi)]j q (d rjjq )— ( -j 
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Jgip ,fi) = exp[-4C A 5(^,/i) + 2A Jg {n h ^)]j g {d mg )— 



V 2 \V 2 J r (%J ' 
where 

r] jq = 2C F Ar{fij,fj) and r) jg = 2C A A T ([i j ,n) . 
The Laplace transforms j(L) are given in Appendix lAl 

4.4 Final resummed distribution in SCET 

With the hard, jet and soft functions in hand, we can now combine them together to form the 
differential inclusive photon distribution. Using the factorization formula, Eqs. ( I54l) - (l55l) with 
the notation of Eq. (Tl5|) . the partonic cross section in SCET takes the form of Eq. (TT9l) : 



d 2 a 
dwdv 



wa(v) H(p T , v, yu) J dkJ(m 2 x — (2Ej)k, ^)S(k, rii, /i) . (90) 



We will now perform this convolution. 

First, note the soft function can depend, in general, on any dimensionless ratio of dot 
products of the directions n^. However, in Section 14.21 we saw that in the formula for the 
resummed soft functions, Eq. ( ITgj) . the only combination which appeared is n 2 ^ 2J ■ This ratio 
can be written in the suggestive form 



2n 12 {2E,jfs _ (2Ej) 2 ^ 
nijn 2 j tu Pt 



Since the soft function is, to all orders, f times a function of k x 2ni2 , if we rescale k — > -£h 

in Eq. ( )90|) . all of the Ej and n« dependence completely disappears, as it must. 

The exact solutions of the RG equations for the soft and jet functions given in Eqs. ( 1781) 
and (189]) involve derivatives acting on a kernel which is just the relevant scale raised to a 
power. Because of this simple form, the convolution in Eq. ( 190]) can be performed analytically. 
In the annihilation channel, the result is 

= wa qq (v) exp [(4CV + 2C A ) S(ji h , fi) - 4C A 5(/i„ /i) + (2C A - AC F ) S(n„ /i)] 

x exp [-2A H (ji ht n) + 2A Ja {n j ,n) + 2A Sm (ji„ /i)] (92) 

/ 2 \ -{2C F +C A )A v (p h ^) , 2 \ (4,C F -2C A )A T (tis,n) 



2 1 fm 2 x \ Vq9 fi-7B%3 



x 17T^\ H ^ PT > u ' Wh^wM 8 ^ — — + d Vm , l i a )— T 
a s {fi) Vt^s m x \ 



>i'i 



where 

Vgg = V J g + Vgg = 2C A A v {n h fi) + (AC F - 2C A )A r (fi s , p) . 
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For the Compton channel 

d 2 a 



dwdv 



wfr qg (v) exp[(4C F + 2C A ) S(ji h , fi) - AC F S(fij, fi) - 2C A S(/i s , fi)] 



x exp [-2A H (fj, h , fj) + 2A Jq {nj,n) + 2A Sqg {fi s , fi)] (93) 

/_2\ -(2C F +C A )A r (ti h ,fi) / 2 



where 



%s = ^ + = 2C F Ar(^-, //) + 2C , j4 A r ( / u s , /i) . 

With these closed form expressions, it is straightforward to evaluate the differential cross 
section numerically. 



5 Cross-checks 

Next, we perform some non-trivial cross checks on the factorization theorem. First, we show 
that in the threshold limit, the expression is RG invariant, that is, independent of fi. Then 
we will show that the singular parts of the fixed order expansion agree with the exact NLO 
differential distribution. We also generate all the plus distribution terms in 1 — w to NNLO, 
and compare to previous results. 

In traditional approaches, one has access to only the renormalization scale fiR, where a s 
is evaluated, and the factorization scale, pif, where the PDFs are evaluated. In the SCET 
approach, there are four scales: the hard matching scale fih, where the hard modes of QCD 
are integrated out, the jet scale, /ij, where the collinear modes are integrated out, the soft 
scale /j, s , where the soft modes are integrated out, and the factorization scale fif, where the 
PDFs are evaluated. In contrast to the scale /xr, the scales coming from the effective field 
theory calculation are all guaranteed to have natural values, since the calculation has been 
factorized into a series of single-scale problems. When combining the RG evolved hard, jet, 
and soft functions with the PDFs, one evolves all of them to a common reference scale fx. 

To check RG invariance, we evaluate the factorization theorem for a fixed scale fih = 
fij = fi s = fif = fi and then show that the cross section is /x-independent. To do so, we 
must expand around the physical, observable machine threshold, M\ — > 0, not the partonic 
threshold mx —> 0, since mx gets integrated over. Recalling Eq. (JHJ) , the invariant mass of 
the hadronic final state is 

M 2 X = m 2 x + &[(l- Xl )v + (1 - x 2 )v] +■■■ 
vv 

near the threshold. Since the SCET operators do not mix, the contributions of different 
partonic channels are separately RG invariant. For the annihilation channel, we have 
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(2E J )fc + ^(l-x 1 ) + ^(l-x 2 



• (94) 



We write the right-hand side in terms of v, which is implicitly a function of y,x%,x 2 and mx 
to make contact with the partonic expressions in the previous section. In this form, we can 
read off that M\ — > enforces m 2 — >■ 0, k — > 0, x\ — > 1 and x 2 — )■ 1 so all the various objects 
approach singular limits. 

Since all the objects are convoluted together in this simplified factorization theorem, it 
makes sense to check RG invariance in Laplace space, which turns the convolution into a 
product. We define the Laplace transformed cross section as 



d 2 a 
dQ 2 dy 



dM x exp 



Ml 



d 2 a 



Q 2 e^ J dM 2 x dy ' 



(95) 



Absorbing e 7E into this definition avoids a proliferation of 7^'s in the Laplace transformed 
expressions. The Laplace transforms of the soft and jet functions are defined analogously (see 
Eqs. (jI2g]l and ([133]) in Appendix [£} . 

For the RG evolution of the PDFs, we can use simplified Altarelli-Parisi equations near 
the endpoint. For the quark PDF, 



dln/i 

The Laplace transform 



2l fq f q /N(x,fJ>) + 2CV7 cusp (a) 



x 1 — X 



(96) 



then satisfies 



f q /N(r, n) 

dfq/N(T,fJ.) 

dln/i 



dxexp 



1 — x 

Te ~lE 



f q/N (x,fj.) 



[2C F 7 cusp In (r) + 2 T / *] f q/N (r, /i) . 



(97) 
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The gluon PDF equation is the same with Cp — > Ca and ^ q —> 7^ 9 . The Laplace transforms 
of the gluon jet function satisfies 



dln/i 



-2CU7CUSP In 



i> 2 



2Y 9 



(99) 



The quark jet function is the same with Ca replaced by Cf and r y J<3 replaced by 7 J<J . The qq 
soft function satisfies 
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din// 99 
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where n. 



'j 



before. 



The other soft function, 

So 



S qg , is the same with 



-4CV + 2C, 



replaced by — 2CU and 7*99 replaced by 7*99. Finally, the hard function satisfies, 



dln/i 



(2C F + cub 



sp In ^ - 7cusp \h(v 2Cf v 2Cf ) + 2 7 
A* 



Hqq(p T ,V,n). (101] 



The function H qg (pT, v, /x) is related to H qq (pT, v, /x) by crossing. Its anomalous dimension can 
be obtained from the above equation by replacing y 2C,F — >■ t) 2CA . 
Putting everything together, RG invariance requires 



dln/i 



H qq -(p T ,v,n)j g (Q 2 ,fi)s 



'/</ 



2£j 



^ fq/N! 



Q 2 - 

Pt 



fq/N 2 



Q 2 



0. (102) 



This equation imposes several constraints on the terms proportional to the universal cusp 
anomalous dimension 7 cttsp , since the /x, p?, Ej, and w-dependence must all vanish. Using 



the above equations and the relationship between rtj,-, Ej and pt in Eq. (19 W . we find that 



For the non-cusp pieces of the anomalous dimension, Eq. f 1 1 2 j) 

(103) 



these constraints all hold, 
implies 

27^9 — 7^ — 7' 5 9? _|_ ry H — Q 

The corresponding relation in the Compton channel reads 

yfi + jfn - Y q - 7 5 " + 7^ = . (104) 

We have calculated all these anomalous dimensions to order a s , verifying these relations. 

Another way to check RG invariance is order-by-order in a s (cf. Ref. [26j for a similar 
example). First by expanding the resummed hard, jet and soft functions, from Eqs. (j6ip . (1891) . 
and (1781) it is easy to check that the matching scales ^hi^j an d fi s cancel between the fixed 
order expansions, such as H(pt,v, fih) and the evolution kernels, such as the eaqp(S(fJ,h, A*)) 
factors in Eq. (|6T|) . To combine the ingredients together to check the overall \i dependence 
at order a s , we would also need a perturbative expression for the PDFs, which is impossible 
since they are non-perturbative. However, since RG invariance should hold with any PDFs, 
a simple trick is to use toy-model PDFs with a convenient form. For example, we can define 
the quark PDF so that f q /^(x, fi ) = 5(1 — x) to all orders at the scale /i . Then, to order a s , 
this PDF at the scale /x is 



fq(x,fj) =5(l-x) + 



An 



-3C F ln 



t4 
» 2 



6(1 -x) + 



in 



-AC, 



■In 4 



1 



x 



(105) 



The gluon PDF is the same with 3Cf — > (3$ and 4C^ — > AC a- Then one can convolute these 
simple PDFs together with the NLO hard, jet, and soft functions to verify /x independence to 
order a s . 

Finally, the fixed order expressions for hard, jet, and soft functions when combined at 
the scales /x^ = fij = fi s = Pt, should reproduce all the singular terms in the exact parton 
level NLO amplitudes from the full standard model. The results for these singular terms are 
presented in Appendix [B] and agree precisely with [3]. In addition, working to order a 2 , we can 
derive all of the terms singular in 1 — w at NNLO. Previous results with NLL resummation [9] 
were able to predict only some of these singularities. 
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Figure 3: Size of the hard and jet function one-loop corrections as a function of the scale for 
different values of px at Eqm =1960 GeV. The right panel shows the optimal scale choice at 
different pr, with the dashed lines denoting our default choice, Eq. (I107p . 



6 Scale choices and matching 

While the resummed result is formally independent of the scales fih, Hj, and fi s , there is 
residual higher-order dependence on these scales if the perturbative expansions of the hard, 
jet and soft functions are truncated at a finite order. To get a well behaved expansion, we want 
to evaluate each contribution at its natural scale, where it does not involve large perturbative 
logarithms. In a fixed order calculation, the presence of several scales can preclude such a 
choice, but since the hard jet and soft functions each only depend on a single scale, we are 
guaranteed that there are scale choices for which large logarithms are absent. 

By examining the form of the resummed distribution, Eqs. (1921) and (193!) . it can be seen 
that the hard, jet and soft scales appear in the cross section only through the combinations 



p\ m\ m\ 



(106) 



Picking fih = pt, fij = mx and fi s = m 2 x /pT as the canonical scales would guarantee the 
absence of large logarithms, but this choice is problematic. To see the problem, recall that 

2 

m\ = —-¥-(1 — w), and the parton-level distribution is singular at w = 1. This singularity 
is integrated over since the hadronic final states are integrated over, and the final photon p? 
spectrum is completely regular. Near w ~ 1, the mass of the partonic final state mx becomes 
small and with the choice fij = mx the coupling constants a s (fij) and a s (/j, s ) are evaluated 
at arbitrarily low scales. Because of the Landau pole singularity of the running coupling the 
convolution integrals are then no longer well-defined. The w ~ 1 part of the integrand is 
suppressed by the resummation, and the contribution from this region of the integral should 
only amount to a power-suppressed correction to the overall result. However, the spurious 
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Figure 4: Determination of \iy On the left is the relative cross section for variations of p,j 
around /ij = pr for _Ecm=1960 GeV. The other scales are chosen to be ph = \if = Pt and 
p s = /A?/ 'fih- On the right, the values of pj which minimize the scale variation at various px's 
are shown for the Tevatron and the LHC. The solid lines show a linear regression to the points, 
and the dashed line is our default choice, Eq. ( I107p . 



power corrections arising in the integration can be of a lower order (and thus of larger size) 
than the physical power corrections to the factorization theorem [58] . 

In [59] it was argued that these spurious singularities are particularly strong in momentum 
space and that it is therefore preferable to perform resummation in moment space. However, 
the effective theory framework allows us to completely avoid the need to evaluate the coupling 
at unphysically small scales. It is not necessary to eliminate the logarithms in the partonic 
cross section, what matters is that the final physical cross section is free of large logarithms. 
Instead of choosing the jet scale pj at the integrand level we should choose the scale after the 
convolution with the PDFs. That is, instead of setting pj = mx, the appropriate jet scale is 
something like the average mass of a jet contributing to the cross section. 

To get a sense of what an appropriate average jet scale should be, let us consider some 
limits. At very large Pt, the relevant scale in the physical cross section is the mass of the 
hadronic final state, so the choice p? ~ M\ = Eq M (1 —Pt/Pt^) i s appropriate. However, at 
moderate px, which is relevant in practice, the appropriate scale choice is less clear. In this 
case, the partonic mass mx at a given p? value can vary kinematically over a large range, 
< mx < Mx, but the fall-off of the PDFs near x — > 1 suppresses the region of large Mx 
and hence of large mx as well. Consequently, the partonic threshold region of small mx 
is enhanced. This dynamical enhancement of was pointed out by [HJ [Tj and was studied in 
detail [20] for the case of Drell-Yan production. It was found that this enhancement is mostly 
effective for relatively high Drell-Yan masses, which corresponds to high px in our case. 

Since we cannot perform the convolution integrals analytically, we will determine the ap- 
propriate choice of pj numerically, following two different procedures. On the one hand, we 
can study the size of the corrections which arise at the different scales. Once the scale is 
chosen appropriately, no large logarithms and associated large corrections should arise. To 
study the size of the corrections, we take the factorized cross section, Eq. (fT9l . as a function 
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Figure 5: Reduction of the factorization scale dependence through matching. The dotted lines 
show the fif scale uncertainty of the unmatched NNLL result, the red lines show the NLO 
uncertainty, and the green band shows the fif uncertainty on NNLL matched to NLO. This is 
for pp collisions at E CM = 1960 GeV integrated over —0.9 < y < 0.9. 

of /i, integrate over the partonic phase space, and compare the tree-level value to the result 
obtained after including the one-loop corrections to either the hard, jet, or soft function. The 
result is shown in Figure The figure shows that the hard corrections are moderate if they 
are evaluated at fi^ ~ px, as expected. The jet function corrections are small at a lower value. 
Looking at the middle panel, we find that the choice fij ~ ^ is reasonable for small pt- For 
larger values of pt, the the optimal scale fij is lower than To be concrete, let us define the 
optimal scale as the scale which minimizes (or in the case of the hard function maximizes) the 
correction. The right-hand panel shows that the choices 



[i-h = Pt , 




(107) 



provide a good approximation to the optimal scale choice as a function of px- For the soft 
scale, we choose fi s = fi^/fih as our default choice and we have checked that the corrections 
are moderate for this scale choice. The plots in Figure [3] are for the Tevatron case, but we 
have also checked that the above scale choices are also valid at the LHC, and that the optimal 
scales for the qq and qg channels are compatible. 

The reasoning behind the above procedure for choosing the scale is that there are no large 
logarithms and thus no large corrections if the scale is chosen appropriately. Another criterion 
for a good scale choice is that the residual scale dependence should be small. To explore this, 
we set fih = A*/ = Pt and fi s = fi^/fih so that the cross section only depends on the single 
scale Hj. We then choose fij such that the distribution is minimally sensitive to variations 
in fij away from its canonical value. In the first panel of Figure H] we show the photon pt 
spectrum integrated over \y\ < 1 at the Tevatron for various values of fij. For simplicity, 
we normalize to the cross section at fij = pt, but since we are only interested in the scale 
dependence, the normalization is arbitrary. The position of the maxima fit nicely along the 
curve fij = 0.56(pt — 1-6^-)- The same procedure at the LHC (14 TeV) gives a best fit 
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Figure 6: Scale variations at the LHC (14 TeV). The lighter bands are NLL and the darker 
bands are NNLL matched to NLO. The unmatched NNLL curves are shown as dotted lines. 



fij = 0.57(pt — l.9-jp—). In the right panel, we show these points, the fits, and our simple 
scale choice, Eq. (110711 . It is comforting that also this criterion leads to similar results. 

So that the results from SCET agree with the NLO partonic cross section in the appropriate 
limit, power corrections must be added through matching. Because of the peculiar kinematics 
of the threshold limit, this must be done with some care. The factorization theorem in SCET 
is derived in the limit where the momentum fractions x\ and X2 of the incoming partons, and 
the partonic threshold variable w, are all close to 1. The resummed cross section is therefore 
only formally /!/ independent for very large px, in contrast to the fixed-order cross section, 
which has additional terms to cancel the \x$ dependence exactly, but only works to order a s . 
These additional terms are not singular in the threshold variables and therefore not reproduced 
by the leading-power factorization theorem. In the phenomenologically relevant regime, x±, x<i 
and w may not be close to 1, and the residual scale dependence might not be small. This NLO 
part of the fif sensitivity can be removed as we match to the NLO partonic cross section, if 
the factorization scale in the NLO cross section is varied appropriately. For the matching, we 
use 



matched / ^2 \ NNLL / ^2 \ NNLL / A% \ NLO 

(108) 



(-,9 \ matcnca / ,9 \ ininliLi / ,o \ injnijIj / i9 \ 

d z cr \ _ ( d z cr \ ( dV \ / dV \ 

dt>dw / V d-ud-u; / V dt>du> J , , , , \ dvdw J 

The subscripts of the last two terms mean set all scales equal to the relevant value of fif. Having 
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Figure 7: Direct photon distributions at the Tevatron, compared to SCET. Green bands are 
scale uncertainty. On the left, comparison is made to CDF data. On the right, the rapidity 
distribution is shown for p T = 200 GeV. The SCET prediction, matched to NLO, is compared 
to the scale uncertainty on the NLO prediction (solid red lines) and to the PDF uncertainty 
(dashed blue lines). 

/if in the matching terms vary in this way significantly reduces the overall fif dependence, 
as can be seen in Figure El This figure also shows that the factorization scale uncertainty at 
large pr is smaller than the uncertainty on the NLO cross section, even without matching. 

With the canonical scales and matching procedure established, we estimate the higher 
order uncertainty by varying the scales by a factor of \ to 2 around their default values. The 
resulting uncertainties are shown in Figure [61 The overall uncertainty is dominated by the 
factorization scale variation. The small bands from variations of fij and fi s should be taken 
with a grain of salt. The above discussion shows that our scale choice is close to the point 
with minimal scale sensitivity, so that the scale variation might underestimate the higher order 
corrections. Also, we observe that the one-loop corrections to the soft function happen to be 
small in our case, much smaller than what was found in other applications. 

7 Results 

To compare to data, we need to deal with the important experimental issue of photon isolation. 
To account for isolation we use the Monte Carlo program jetphox. This program includes 
both the NLO partonic cross section and a fragmentation contribution, applying a user-defined 
isolation criteria. To correct the SCET distributions for isolation, fragmentation, and finite 
NLO effects, we match to jetphox, i.e. we use the output of this program for the NLO cross 
section in the matching relation Eq. (JTUHJ) . To compare to the DO data [60], we attempt to 
match their isolation criterion by demanding less than 10% of the energy in a cone of R = 0.4 
around the photon be hadronic. For the CDF data [BTl E2], we require less than 2 GeV of 
energy inside the R — 0.4 cone. Some studies of sensitivity to isolation parameters can be 
found in [6TJ and we do not attempt to reproduce them here. 

In addition, we apply to all the Tevatron theoretical calculations an overall rescaling 
of 0.913 (taken from [6H [62]) to account for underlying event, multiple interactions, and 
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Figure 8: Fixed order and resummed comparison to DO and CDF data. Left plots show 
the LO and NLO scale uncertainties. Right plots show the SCET predictions with NLL 
resummation or with NNLL resummation matched to fixed order. The dashed blue lines are 
PDF uncertainties. The curves are all corrected for isolation, fragmentaion, and hadronization 
as described in the text, while the reference distribution dcr^o is the fully inclusive NLO 
distribution without corrections. 



hadronization. The DO data corresponds to 380 pb _1 of integrated luminosity at -Ecm = I960 
GeV, integrated over —0.9 < y < 0.9. The CDF data corresponds to 2.5 fb _1 of integrated 
luminosity at Eqm = I960 GeV, integrated over —1 < y < 1. For all calculations, including 
jetphox and scale uncertainties, we use the MSTW 2008 NNLO PDFs [63]. The rationale 
behind this choice is that our calculation includes the dominant NNLO corrections. 

The scale uncertainties for the fixed order result include variation of the factorization 
scale Hf, the renormalization scale and a fragmentation scale M' . The fragmentation 
scale is related to collinear singularities in final state photon emission from, for example, qq 
final states, which are relevant starting at NLO. For simplicity, we call all these scales \i 
and vary them together. We define the NLO uncertainty as the maximum and minimum 
value of the prediction from varying these scales between |p r < /i < 2pr- For the SCET 
prediction, we vary the jet, hard, soft and factorization scales. The largest uncertainty is 
from the factorization scale variation, even after the proper matching to NLO (see previous 
section), and so we use the fif dependence for the SCET uncertainty bands. Again, we take 
the maximal and minimal values along the range ^p? < fif < Ipt- 
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Figure 9: Predictions for the inclusive direct photon distribution at the LHC. Left plots 
show the LO and NLO scale uncertainty. Right plots show the SCET predictions with NLL 
resummation or with NNLL resummation matched to fixed order. The dashed blue lines are 
PDF uncertainties. No correction for isolation or hadronization is included. In contrast to 
Figure EJ here NLO refers to the inclusive direct photon distribution whose central value is 

(dir) 
NLO' 



identical to the reference distribution d<rf dl 



Figure [7] shows the px and rapidity distributions at the Tevatron. The px distribution is 
compared to CDF data [61] [62j and the rapidity distribution only to the inclusive NLO result 
and the PDF uncertainties. No comparison to data has been made in the rapidity plot because 
all of the published Tevatron data differential in the photon rapidity is differential in the jet 
rapidity as well, for which our factorization theorem does not apply. Nevertheless, such a 
comparison would be interesting as there is a significant discrepancy between the SCET result 
and the NLO prediction. 

For more detail, we show in Figure [S] the normalized pt spectra and compare to CDF |JU 
16*2"] and DO data [HO]- In this figure and in the LHC plots in Figure |H1 we normalize to 
<7 NLO , the inclusive NLO direct photon cross section, without isolation cuts and fragmentation 
contributions, evaluated with the default scale choices. The left plots show the LO and 
NLO distributions, matched to jetphox, with the blue dashed lines indicating NLO PDF 
uncertainties (from the MSTW 2008 NNLO PDFs). The right plots show the predictions from 
SCET at NLL and NNLL, also matched to jetphox, with the appropriate PDF uncertainties 
included as well. Note that at high px, the scale uncertainty for the SCET result is smaller 
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than the PDF uncertainty, while for the jetphox, NLO result, it is not. 

Figure 13 shows the SCET prediction at the LHC, with E CM = 14 TeV and E C m = 7 TeV, 
integrated over — 1 < y < 1. The two energies give results that are qualitatively very similar, 
but of course, the higher energy machine would be capable of producing higher px photons. 
Note the large PDF uncertainties at high px- These distributions indicate that the quark, 
antiquark, and gluon PDFs at large x can be measured effectively with direct photon produc- 
tion at the LHC. The PDF uncertainties are slightly larger with the SCET cross section than 
with the NLO cross section, due to slightly different scales and x-values where the PDFs are 
evaluated. 

8 Conclusions 

We have shown how to resum the direct photon distribution using Soft-Collinear Effective The- 
ory. This is the first physical process calculated in SCET involving more than two collinear 
directions. The factorization theorem we derived passes a number of non-trivial checks. In 
particular, renormalization scale independence arises only after a cancellation of the depen- 
dence on angles rii-rij appearing in the soft function against s, t, and u dependence in the hard 
function. We have calculated all the relevant objects to one loop, showing that the dependence 
on the various kinematic variables is indeed of the right form for the factorization theorem to 
hold. 

Our closed-form expression for the distribution allows for straightforward numerical inte- 
gration and comparison to data. The agreement with Tevatron data is very good, although 
the comparison is complicated by the issue of photon isolation. However, the factorization 
theorem holds at large pt, where the isolation is less relevant. We also show results at higher 
center-of-mass energy, relevant for the LHC. There, we have found that the theoretical un- 
certainty is much smaller than the PDF uncertainty. Thus, the resummed direct photon 
distribution will be an effective tool for measuring the PDFs at the LHC. In addition, there 
is a significant difference between the SCET prediction at high pt at the LHC and the NLO 
prediction. This is not surprising as there are large logarithms in this region which SCET 
resums to all orders. In particular, our NNLL resummed result has all the singular parts of 
the partonic cross section at NNLO (except for the ^-function part, which is only incompletely 
known). Based on experience with other processes, such as Drell-Yan or Higgs production, 
our result is expected to provide a good approximation to the full NNLO cross section. 

Besides being of phenomenological importance, the calculations in this paper are easily 
generalizable to other fundamental processes at hadron colliders. The obvious example is W 
or Z production at high pt- For this case, the factorization theorem is identical. The jet 
and soft functions are also the same, and the hard function is the same up to corrections 
finite in mz and my/- Since W/Z production is free of the complication of photon isolation, 
it is cleaner phenomenologically. We are currently pursuing the analysis and will present 
our results elsewhere. For direct photon production, it would also be interesting to treat 
the photon isolation cut in the effective theory, since it is known that the corresponding 
perturbative expression involves large logarithms which make the fixed order calculation of 



33 



this effect problematic [I]. 

The next steps in complexity are to consider distributions which are hadronically no longer 
fully inclusive, and to include outgoing jets for the proton remnants. In either case, a modified 
factorization theorem is needed. Including proton remnant jets is necessary to get away from 
the end-point region where the leading partons carry almost all of the proton momentum. 
To get to the phenomenologically more interesting region of moderate momentum fractions, 
one needs to account for the energetic partons down the beam pipe, a problem which has 
only recently been considered in SCET [51] . Given the progress in this field over the past few 
years, it now becomes possible to analyze also complicated collider processes with effective field 
theory, in particular processes with several hadronic jets and nontrivial kinematical restrictions 
on the hadronic final state. 
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A Fixed Order Expansions 



In our notation, all the anomalous dimensions are expanded as series in 
The QCD (3 function is 
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The four- loop coefficient f3$ is known as well (65j [66], and can be found for example in 
The cusp anomalous dimensions is 
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The anomalous dimensions describing the evolution of the quark and gluon PDFs near 
x = 1 are 
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The three loop splitting functions were calculated in [671 EE] . Explicit expressions for anoma- 
lous dimensions 7/ and 7/ at three loops can be found in [20j [23]. 
The hard function can be written as 



Pt 



H(pr,n,v) = h(\n—,v) . 



To order a 2 it is 
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As described in Section HJ the anomalous dimensions can be extracted from the general result 
J. Explicity, 
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The three loop anomalous dimension for the hard function is just 7^ = 272 + 7I , where the 
anomalous dimensions 7 9 and 7 9 were defined and given to three loops in [37]. 

The u-dependence is different in the two channels. It shows up in the functions and 
the constants c¥ . For the annihilation channel, we find 
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and for the Compton channel, 
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The solution of the RG equation for the jet function is given in terms of its Laplace 
transform 
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This Laplace transform has a perturbative expansion in a s which only depends on the dimen- 
sionless ration Q 2 //i 2 . Writing i(Q 2 ,/i) = j(L), where L = log%, the expansion becomes 
very similar to that of the hard function 
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For the quark jet 
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The three loop anomalous dimension ^ an d the two loop constant c J 2 q are also known [55], 
but are not necessary for NNLL resummation. For the gluon jet 
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j ( 331153 6217vr 2 _ 583tt 4 647r 2 C 3 \ . . 



1096 llvr 2 \ „ „ /368 4vr 2 



/ 42557 2612 I6C3 , 154vr 4 \ 2 /3622 80tt 2 448C 3 \ ^ 2 ^ 2 
+ {^m " "243" " ~W + "l35~ J C ^ /Tf + ^729" + ~81 2^T J 

+ J4145 _ 4^ _ 60^ _ ^ _ 2C%nfTp _ U CF<n 

The soft function S(k,fi,riij) depends in addition to the scales k and fi on the angles 
riij = rii ■ rij between the Eikonal lines. However, this dependence must be universal for the 
factorization theorem to hold. The Laplace transformed soft function 

s(K,[i,n ij ) = J dke^ f--^-J S(k,fi,nij) (133) 

has a perturbative expansion in a s which only depends on one dimensionless ratio. We can 
write s(/t, /i,n,j) = s(L) where 
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The expansion is now similar to the hard or jet functions 
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The coefficients in the annihilation channel are 
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The Compton channel is identical, but with the prefactor Cf — \Ca replaced by \Ca- This is 
a consequence of Casimir scaling, which holds at least to three-loop order. 

B NLO and NNLO singular terms 

To obtain the singular terms, the resummed results, Eqs. (|9"2|) and (1931) . should be expanded 
order-by-order in a s . To do so, we set all the scales equal fih — fij — // s = fif — In the 
limit of equal scales, the various evolution factors S(v,fi) and A 7 (i/, //) and the quantities r)q q 
and r) qg all vanish. Before setting the scales equal, we expand the kernel using 
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and perform the derivatives with respect to rj. Then we take rj — > 0. The resulting expressions 
are lengthy. To save space we only quote the result for fi = p^, for which we find 
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This formula holds for either channel, with the appropriate cr(v) and hard, jet and soft function 
coefficients. The coefficient A2 is not completely known, so we do not include our partial 
results. We have checked that the 0(a s ) results agree with j3] and the a 2 s ln -/l""^ and 

results agree with [9]. We find a small discrepancy with these authors for 

term in the Compton channel, but otherwise we confirm their results. The 
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